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ABSTRACT 


The NASA Dryden Flight Research Center is exploring the optimum thrust-vector angle on aircraft. 
Simple aerodynamic performance models for various phases of aircraft flight are developed and 
optimization equations and algorithms are presented in this report. Results of optimal angles of thrust 
vectors and associated benefits for various flight regimes of aircraft (takeoff, climb, cruise, descent, final 
approach, and landing) are given. Results for a typical wide-body transport aircraft are also given. The 
benefits accruable for this class of aircraft are small, but the technique can be applied to other 
conventionally configured aircraft. The lower L/D aerodynamic characteristics of fighters generally 
would produce larger benefits than those produced for transport aircraft. 
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quadratic equation coefficients 

acceleration along the flightpath axis 

drag coefficient 

minimum C D 

lift coefficient; W/(qS) 

lift coefficient caused by flap deflection 

C L at C D 0 
drag 

aerodynamic efficiency, lift-to-drag ratio, E = L/D = c l /c d 

rolling or braking friction force 

flightpath angle 

forces along the flightpath axis 

forces perpendicular to the flightpath axis 

acceleration caused by gravity 

efficiency factor in quadratic drag equation 

lift force 

lift-to-drag ratio 

lift force at rotation 

dynamic pressure 

range (takeoff or landing distance) 

area of wing, 3500 ft 2 

time 

thrust 
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touchdown 
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^rotation 


to be determined 
aircraft speed 
calibrated airspeed 
aircraft speed at touchdown 
aircraft weight 
weight on wheels 
flightpath axis 

perpendicular to the flightpath axis 
flightpath angle 

angle between the flightpath angle and the thrust vector 
pitch attitude at rotation 

coefficient of rolling friction (can include braking effect) 


INTRODUCTION 


One task an aircraft designer always has in the design of conventional aircraft is determining the 
inclination of the thrust angle to the fuselage reference line. The optimal inclination is a function of 
aircraft characteristics and flight phases; and for aircraft with no variable thrust-vectoring capability, the 
selected inclination angle must be a weighted compromise considering all relevant factors. 

The concept behind the optimal thrust inclination angle is simple: if the angle is positive (thrust is 
tilted down), an additional vertical force is contributed that decreases the lift requirements of the wing. 
However, this inclination also slightly reduces the horizontal thrust component; thus, an optimization 
procedure (refs. 1 and 2) is required to determine the optimal angle. This report investigates the 
relationship of thrust inclination angle for various flight phases. The results can provide insight to the 
designer for selecting the thrust-vector angle. 

Simplified equations of motion are also included for the various flight phases. These sets of equations 
are then optimized with respect to the thrust inclination angle; the benefits of optimization also are 
derived from the analysis. The flight regimes investigated are: 

• takeoff (minimum takeoff distance). 

• climb (maximum rate of climb for a fixed thrust level). 

• cruise (minimum thrust; minimum fuel consumption). 

• descent (best glide range for a fixed thrust level). 

• final approach (minimum approach velocity). 

• landing (minimum rollout distance). 
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The optimization is performed for the full quadratic drag representation, and comparisons are made 
with a linear drag representation. In addition, an analytical optimization procedure is developed for 
steady-state flight conditions and the landing rollout. Example results are given for a representative wide- 
body transport aircraft. 

THRUST ANGLE OPTIMIZATION 

The ideal inclination of the thrust angle is a function of the particular flight phase. For aircraft with no 
thrust-vectoring capability, the selected fixed inclination angle depends on the aircraft type and mission; 
and the selected angle is weighted towards the dominant flight phase. For the few aircraft that currently 
have— and the future aircraft that will have— thrust-vectoring capability, optimization will change as 
various flight phases and conditions are encountered. 

The equations of motion (ref. 3) are intricate, especially because of the quadratic representation of 
drag required to assure accurate results. In addition, takeoff and landing conditions are dynamic and thus 
require time-domain solutions; whereas steady-state flight phases can be analyzed with algebraic 
solutions. The equations and analysis developed apply equally well to conventional aircraft of all sizes, 
however, actual results are affected by the aerodynamics of the particular aircraft. Lift and drag 
characteristics are the most dominant aerodynamics affecting results. The results of the techniques 
developed herein will be illustrated with a wide-body transport for the different flight phases. The 
analysis identifies the optimal thrust-vector angle and the accruable benefits. Table 1 shows the flight 
conditions and required aerodynamics for the flight phases analyzed, which are representative of a typical 
wide-body transport. Analytical optimization procedures are developed in appendices A and B for 
steady-state flight conditions and the landing rollout. 

OPTIMAL THRUST ANGLE AND BENEFITS FOR TAKEOFF 

The optimization objective for the takeoff flight phase is to minimize the distance to aircraft rotation. 
Figure 1 shows the forces acting on an aircraft at the center of gravity during takeoff acceleration. The 
resulting equations of motion for takeoff with a thrust-vectoring capability, r \ , can be written as follows: 

^ F x = 0 = -D + 7cos(r|) - WA x /g - F (1) 
^F z = 0 = L + rsin(ri)- W + WoW . (2) 



Figure 1. Aircraft forces during takeoff acceleration to rotation. 
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Table 1 . Aircraft parameters and flight conditions. 
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Takeoff 

Climb 
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Final 
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460,000 
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280,000 

280,000 

280,000 

T 

100,000 

TBD 

TBD 

TBD 

TBD 

-70,000 

Cd o 

0.050 

0.018 

0.018 

0.018 


0.100 
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0.105 

0.105 

0.105 



C K 


0.16 

0.16 

0.16 



s 

0.3 






Cr 

1.300 at 

0.447 

0.459 

0.228 

1.000 

0.000 


rotation 






C D 


0.0267 

0.0274 

0.0185 

0.1200 

0.100 

E 


16.80 

16.80 

12.30 

8.33 


V 

0.02 





0.40 

0 , 

10.00 






rotation 








Solving for the weight on wheels, WoW , from equation (2) yields: 

WoW = W-L-T sin(ri). (3) 

During the horizontal portion of the takeoff acceleration, the lift, L, primarily is caused by flap deflection. 
The wheel/runway friction force can then be expressed as: 

F f = \JlWoW . (4) 

Rewriting equation (1) using equation (4), 

A x = g|Tcos(r|) - D - \\.WoW]/W , (5) 


4 



where the drag variation up to the point of rotation is 


D = qS 




( 6 ) 


The velocity at which the aircraft can take off is the L rotat i on criterion defined by 
^rotation ^-TsinOl-h 0 rotation ) , where L rotation is the total aircraft lift that would be generated if 
the aircraft is rotated to the desired pitch attitude, 0 rotation , and is a function of velocity defined by 


L rotation ^L’ 


(7) 


where C L is a function of angle of attack at the point of rotation. 

Equations (5) and (6) cannot be analytically solved for the minimum range, R\ time-series solutions 
must be calculated as follows: Selecting a set of aircraft conditions and a nominal runway friction 
coefficient, |X ; A x (eq. (5)) is calculated and integrated twice to obtain V and R (eqs. (8) and (9)) for a 
selected thrust-vector angle, r\: 


A x = g[Tcos(r|) - D - \iWoW]/W (5) 

V = \A X dt (8) 

R = JV dt (9) 

Equations (6), (3), and (7) are calculated to obtained updated D, WoW , and L rotation values for the new 
velocity: 


D = qS\ 


C D +k(C L 


~ C 0 T 


( 6 ) 


WoW = W-L 


rotation 


- T sin(Ti) 


(3) 


rotation 


= qSC 


L ' 


(7) 


The L criterion is also calculated at each instant of time: 

rotation 

^rotation — ^ — 7"sin (T| + ^rotation) ' When the L ro,ation criterio " is raet - ,he airCraft haS rCached 

takeoff velocity, thus defining the takeoff distance. (In practice, rotation could start before the L rotalion 
criterion is met because time is needed to rotate the aircraft and during that rotation period, the aircraft 
has accelerated to a faster speed than required for liftoff.) This analysis is repeated for a range of r| • 
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Figure 2 shows results representative of a wide-body transport for a range of p ; a typical p for a 
concrete runway is 0.02. For a given p , R reduces as the thrust-vector angle increases from 0° to 
approximately 12°, at which point the minimum takeoff distance is achieved. The thrust-vectoring 
benefits are modest; R is reduced approximately 200 ft for p = 0.02 . These results are similar for the 
range of p presented; however, the benefit of thrust vectoring does increase slightly with increasing p . 
The effect of a nominal p of 0.02 increases takeoff roll approximately 700 ft over the p = 0 condition; 
accelerating the aircraft mass with the given level of thrust is the primary factor in determining the 
takeoff roll distance. 



Figure 2. Variation of takeoff distance with thrust- vector angle for various coefficients of rolling friction. 

OPTIMAL THRUST ANGLE AND BENEFITS 
FOR CLIMB, CRUISE, AND DESCENT 

The optimization objectives selected for the climb, cruise, and descent flight phases are the maximum 
rate of climb for a fixed thrust level, minimum thrust (minimum fuel consumption), and the best glide 
range for a fixed thrust level, respectively. 

Because these different flight phases are all steady-state flight, the forces acting on the aircraft can be 
generalized (fig. 3). The resulting equations of motion for these flight conditions with a thrust-vectoring 
capability, r[ , and a flightpath angle (FPA), y , can be written as follows: 

V F x = 0 = -D + Tcos(q) - Wsin(Y) ; (10) 

= 0 = L + T s\n(r]) - W cos(y) ■, (11) 

where 1) is measured relative to the flightpath. 
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Figure 3. Aircraft forces during steady-state flight phases of climb, cruise, and descent. 

Analysis With Quadratic Drag Formulation 


Drag can be represented as 

D = qS\ [ C Do + k(C L - C L o ) 2 ] . (12) 

Rearranging equation (11) gives 

L = [Wcos(y) - rsin(ri)] = qS C L \ (13) 

substituting C L from equation (13) into equation (12) gives 

D = qS^C D +M(Wcos(y)-rsin(TO)/(gS)-C^} 2 J. O 4 ) 

Rearranging equation (10) gives 

D = -Wsin(Y) + 'Tcos(r)). (15) 

Substituting equation (14) into equation (15), 

qs[c D + &{( W cos(y) - T sin(r)))/ (qS) - C^} 2 J = TcosOl) - Wsin(Y). (16) 

Expanding equation (16) and collecting thrust, T, terms in the form 

aT 2 + bT + c = 0 (17) 
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produces the coefficients 


a = sin(ri) 2 ; 

b = -2Wcos(y)sin(ri) + 2qS C , sin(ri) -qS cos (r\)/k\ and 

O 

c = W 2 cos(y) 2 - 2qSWC L cos(y) + q 2 S 2 C L 2 + q 2 S 2 C D ^/k + qSWs\n(y)/k . 

Selecting a set of aircraft conditions including FPA, the quadratic equation (17) is solved for thrust, T, 
for specific thrust-vector angles, T). Only one solution of the quadratic equation provides meaningful 
results. (For the examples calculated, the thrust produced by the other root is of the order of 10 at 
■q = 0.) Plotting thrust as a function of thrust-vector angle will provide the optimal thrust-vector angle 
for minimum thrust required for steady-state flight at the selected FPA. These results apply to climb, 
cruise, and descent flight conditions, depending on the FPA selected. 

Climb Conditions 

Figure 4 shows results for a nominal climb FPA of 3°. For this nominal FPA, a minimum thrust 
condition is achieved at a thrust-vectoring angle of approximately 3.4° and results in a thrust reduction of 
approximately 90 lbf. Also plotted are data for a FPA of 3.01 15° (fig. 4). A comparison of the two sets of 
data shows that for a constant thrust setting, optimal thrust vectoring increases the FPA 0.0115°, 
a 0.38-percent increase. 



Figure 4. Variation of thrust required with thrust-vector and flightpath angles in a climb. 
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Cruise Conditions 


Figure 5 shows results for a level cruise flight condition. For this cruise condition, a minimum thrust 
condition is achieved at a thrust-vectoring angle of approximately 3.6° and results in a thrust reduction of 
approximately 40 lbf. 



Thrust-vector angle, deg 
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Figure 5. Variation of thrust required with thrust-vector angle for cruise flight. 


Although the benefit is small for the cruise flight phase, the takeoff-to-landing benefit is greater, 
percentage-wise, than the instantaneous benefit at any one flight condition. This greater benefit is caused 
by the cumulative effects of reduced fuel consumption for the entire cruise portion of flight and also 
includes the fact that not as much fuel has to be loaded for takeoff. This same cumulative benefit applies 
to the climb flight phase, although the benefit is smaller because of the short amount of time required to 
climb to cruise conditions. The beginning-to-end calculation of benefits is beyond the scope of this 
report; however, a representative penalty incurred carrying unnecessary fuel is approximately 4 Ibm/hr 
for each 100 Ibm of fuel carried. 

Descent Conditions 

Figure 6 shows results for a nominal descent FT A of —3°. For this nominal FPA, a minimum thrust 
condition is achieved at a thrust-vectoring angle of approximately 0.8° and results in a thrust reduction of 
less than 1 lbf. Optimization is minimal at this flight condition because of the fact that the thrust level 
required for this descent flight condition is low. As the thrust level approaches 0 lbf, the benefit of the 
vectoring thrust approaches 0 percent. The best glide range for any aircraft with 0-lbf thrust is 
proportional to the lift-to-drag ratio, L/D , and is defined by the relation y ~ -1 / ( L/D ) . 
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Figure 6. Variation of thrust required with thrust-vector angle in a descent. 


Analysis With Simplified Drag Formulation 

In situations where the aircraft is operating near its best aerodynamic efficiency, the L/D can be 
used for obtaining qualitative results. Representing L/D by E, an expression for drag is D = L t. 
Substituting this expression into equation (10) and solving simultaneously with equation (11) will 
produce thrust as a function of the aircraft weight; FPA; thrust-vector angle; and E: 

T = W[sin(v) + £ _1 cos(y)]/[cos(ti) + E 1 sin(r|)]. (18) 

Figures 7-9 show the results of this simplified formulation compared with the previous results for the 
climb, cruise, and descent results (figs. 4-6). 

For the climb condition (fig. 7), the complete and simplified results are essentially identical. For the 
cruise condition (fig. 8), the simplified result provides a very good approximation and indicates onl> a 
slight reduction in both the optimal thrust-vector angle and the thrust required. For the descent condition 
(fig. 9), the results are significantly different: although the complete analysis indicated essentially no 
benefit for thrust vectoring, the simplified analysis indicates a thrust reduction benefit of 23 Ibf at a 
thrust-vectoring angle of approximately 4.7°. The difference for the descent flight condition is mainly 
caused by the fact that the linear representation of drag as a function of lift is significantly different than 
the quadratic drag representation. The linear representation of drag agrees with the quadratic 
representation of drag only near the maximum L/D for the aircraft, which is generally where aircraft arc 
designed to cruise (thus the good agreement for the cruise flight condition). 


10 




5.1 


FPA =|3.0 
Quadratic Ci 



Figure 7. Comparison of climb thrust required with thrust-vector angle for a quadratic and linear drag 
representation. 
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Figure 8. Comparison of cruise thrust required with thrust-vector angle for a quadratic and linear drag 
representation. 
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Figure 9. Comparison of descent thrust required with thrust-vector angle for a quadratic and linear drag 
representation. 

Analytical Optimization 

For the simplified formulation, the optimal thrust-vector angle can be analytically determined for the 
climb, cruise, and descent flight phases. Appendix A shows the development of the optimal angle. The 

resulting optimization relationship for thrust-vector angle is tan(T|) = 1/E. 

Interestingly, the parameter E totally defines the optimization result. These analytical results are 
identical to the optimal thrust-vectoring angle obtained from the simplified analysis (table 2). The results 
from the complete analysis are similar for the climb and cruise flight phases, but are significantly 
different for the descent phase. 


Table 2. Comparison of optimal thrust-vector angle, 
rj, for the various analysis techniques. 



Flight phase 

Climb, 

deg 

Cruise, 

deg 

Descent, 

deg 

Complete analysis 

3.4 

3.6 

0.8 

Simplified analysis 

3.4 

3.4 

4.6 

Analytical expression 

3.4 

3.4 

4.6 


12 




OPTIMAL THRUST ANGLE AND BENEFITS FOR FINAL APPROACH 


The optimization objective selected for the final approach flight phase is to minimize approach 
velocity, which will result in a reduced rollout distance when touchdown occurs. The steady-state 
equations of motion for the final approach are the same as those in the previous section for climb, cruise, 
and descent. The difference is that in the previous examples, maximizing aerodynamic efficiency was the 
objective; whereas in the final approach, a minimum touchdown speed is the objective. This minimized 
touchdown speed will be obtained at a higher thrust level than would be required for an approach 
condition with no thrust vectoring. 

Repeating equations (10) and (11), ^F x = 0 = -D + Tcos(r[) - W sm(y)- 

F z = 0 = L + T sin(ri) - Wcos(y). Solving each equation for r\, 

cos(ri) = [D + Wsin(y)]/T ; (19) 

sin (r| ) = [Wcos(y) -L]/7\ (20) 

Dividing equation (20) by equation (19) yields: 

tan(r|) = [Wcos(y) - L]/[D + WsinM] - (21) 

The intersections of any two of the three equations ( 1 9)— (2 1 ) for various y and thrust levels will 
provide trimmed flight conditions. However, because equation (21) is not a function of thrust, it provides 
a single curve, a locus of trimmed solutions as a function of y. Equations (21) and (19) are then used to 
obtain thrust-vector angle optimization. 

Figure 10 shows aircraft velocity as a function of thrust-vector angle and thrust level for y = -3.0° ; 
trimmed flight conditions are indicated by the dashed line. In general, increased thrust levels clearly will 
permit significant decreases in aircraft velocity. For the descent conditions with no thrust-vectoring 
angle, the velocity is approximately 279 ft/sec at a thrust level of approximately 18,900 lbf. However, if 
the aircraft were trimmed at a thrust level of 60,000 lbf and a vectoring angle of approximately 78°, the 
aircraft could be flown at a stabilized velocity of approximately 248 ft/sec; more than 30 ft/sec slower 
than without vectoring and a more than 10-percent reduction in approach speed. 

Figure 11 shows that at low vectoring angles, the aircraft can be trimmed at the same thrust 
(18,900 lbf) but with a vector angle of 13.5°, thus reducing approach speed by a little more than 2 ft/sec. 
At a vector angle of approximately 7°, the thrust can be reduced approximately 130 lbf and still maintain 
the same FPA. 
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Figure 10. Variation of approach speed with thrust-vector angle and thrust. 



Velocity, ft/sec 
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Figure 11. Variation of approach speed with thrust-vector angle and thrust. 





OPTIMAL THRUST ANGLE AND BENEFITS FOR LANDING ROLLOUT 


The optimization objective selected for the landing rollout flight phase is to minimize rollout distance 
when the aircraft has touched down. Figure 12 shows the forces acting on an aircraft during the landing 
rollout. The steady-state equations of motion for the landing with a thrust-vectoring capability, r| , can be 
written as follows: 


£F x = 0 = -D + T cos(t)) - WA/g - F f ; (22) 

^F z = 0 = L + Tsin(r\)-W + WoW. (23) 

These equations are identical to the takeoff roll equations (eqs. (1)— (2)). Assuming that the lift, L, is 0 lbf 
because of spoiler and speedbrake deflection, solving for the weight on wheels, WoW , from 
equation (23) is expressed as 


WoW = W-rsin(ri). (24) 

The drag variation during the rollout can be expressed as 

D = qSC n . (25) 

Positive thrust is defined as positive in the forward direction; therefore, for the landing rollout condition 
where thrust will be reversed, thrust ( T ) will have a negative value. This definition implies that when the 
thrust is reversed, a positive r) will provide a downward thrust component, thus increasing the rolling 
friction force as 


F f = p(WW). 


(26) 



Figure 12. Aircraft forces from touchdown and through deceleration to a stopped condition. 


Rewriting equation (22) and solving using a time-series method is calculated as follows: Selecting a 
set of aircraft conditions and a runway surface coefficient, p ; equations (27), (28), (29), (25), and (24) 
are calculated for a selected thrust- vector angle, Tp 

A x = g[Tcos{r()-D-\iWoW]/W\ (27) 
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V - V touchdown + \ A x dt ’ 


(28) 


R = jvdf, (29) 

D = qS C D ; and (25) 

WoW = W-rsin(Ti). (24) 

The landing rollout starts with a touchdown velocity equivalent to a calibrated airspeed of 120 kn. The 
time at which the forward velocity equals 0 kn defines when the aircraft has stopped. This analysis can 

repeated for a range of r\. 

Figure 13 shows rollout distance results as a function of thrust- vector angle and p for a wide-body 
transport at touchdown conditions. The thrust is -70,000 lbf (reversed); = 0.10 ; and the nominal 
p = 0.4 (which includes the effect of near maximum braking capability). For this configuration, the 
optimal thrust-vector angle is approximately 22° and has reduced the rollout distance approximately 26 ft 
(2.7 percent). Also shown is the effect of p , which is the most important factor in decelerating the aircraft. 

Figure 14 shows rollout distance as a function of time to stop for the same data as figure 13 with 
p = 0.4 , and shows that the time to stop is reduced by approximately 0.26 sec (2.8 percent). In addition 
to the benefits of optimization being small, the rollout distance and time to stop are quite short. This 
shortness is because of the model used, which assumes the application of the friction caused by braking, 
the drag effect caused by spoilers, and reverse thrust all are applied at the instant of touchdown. 
Commercial aircraft typically have many seconds after touchdown before these various items are 
activated. In any case, the incremental changes indicated by this analysis are representative of actual 
results. Figure 15 shows the effect of C D and thrust. The C D has a minor effect on rollout distance; 
thrust has a significant effect on rollout distance. 

Appendix B develops the following simplified optimization expression for ri , which produces the 
minimum landing rollout distance: tan (r| ) = p. These results are identical to the optimal thrust- 
vectoring angle obtained from the complete analysis (table 3). 


Table 3. Comparison of optimal thrust- vector angle, 
r| , for various analysis techniques. 
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Figure 13. Variation of rollout distance with thrust-vector angle and coefficient of braking friction. 
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Figure 14. Variation of rollout distance and time to stop with thrust-vector angle. 
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Figure 15. Variation of rollout distance with thrust-vector angle, C D ^, and thrust, T. 


CONCLUDING REMARKS 


Analyses were performed to determine optimal thrust-vector angle as a function of various flight 
conditions. The flight regimes explored were takeoff, climb, cruise, descent, final approach, and landing. 
The equations developed for the various flight phases are simplifications of actual conditions, but the 
resulting relationships accurately depict the optimal thrust-vector trends and provide an indication to the 
aircraft or propulsion controls designer. Results for a typical wide-body transport are as follows: 


Phase 

Optimization 

Improvement, 

percent 

deg 

Takeoff 

Distance 

Reduced 2.5 

12.0 

Climb 

Rate of climb 

Increased 0.4 

3.4 

Cruise 

Thrust 

Reduced 0.2 

3.6 

Descent 

Glide range 

No improvement 

0.8 

Final approach 

Velocity 

Reduced =10.0 

= 78.0 

Landing rollout 

Distance 

Reduced 2.8 

22.0 


Analytical optimization developed two expressions for the optimal thrust-vector angle: one as a 
function of runway friction coefficient for the landing condition, and the other as a function of the 
lift-to-drag ratio, L/D , for all the “up-and-away” flying conditions. These results compared very well for 
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the climb, cruise, and landing rollout flight phases. The process easily can be extended to complex 
aircraft configurations and flight conditions. 

Generalized observations for the various flight phases are as follows: 

• Takeoff: Thrust-to-weight ratio dominates, and increasing values give faster acceleration and 
increased optimal thrust-vectoring angles and associated benefits. 

• Climb, cruise, and descent: Lower L/D configurations benefit most from thrust vectoring 
because of larger optimal thrust-vector angles. 

• Final approach: High thrust-to-weight configurations benefit most, assuming very large thrust- 
vectoring angles can be achieved. 

• Rollout: The runway friction coefficient is dominant, although large levels of reversed thrust 
would be beneficial. 


Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, March 7, 2000 
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APPENDIX A 


OPTIMAL THRUST ANGLE FOR STEADY-STATE FLIGHT 


The steady-state equations of motion for climb, cruise, and descent flight phases with a 
thrust-vectoring capability, q , and a flightpath angle, y, can be written as follows: 

^ F x = 0 = —D + T cos(ri) - W sin(y) ; (A* 1 ) 

'Jj F z = 0 = L + T sin(ri) - Wcos(y) . (A-2) 

Rewriting equation A-2, 

L = Wcos(y) - T sin(ri). (A-3) 

Defining E = L/D ; thus, 

D = L/E. (A-4) 

Substituting expression (A-4) into equation (A-l) and then substituting equation (A3) inio 
equation (A-l) results in 

0 = -WE -1 cos (y) + T£ _1 sin(q) + E cos(ri) - IV sin(y) . 


Solving for T, 

T = fT[sin(y) + £ _1 cos(y)]/[cos(ri) + E 1 sin(r|)] . (A-5» 

To find the minimum thrust, T, with respect to q , the necessary condition is that the partial derivative 
dT/dr\ = 0 ; thus: 


dT/dr\ = 0 

= [(cos(q) + £ _1 sin(q)} {0} - W{sin(y) + £ _1 cos(y)}{-sin(q) + E~ l cos(q) } 1 (A-6t 

_i 2 

/[cos(q) + E sin(q)] . 

Upon reducing and assuming { sin (y ) + E cos ( y ) } ^ 0 , 

0 = -sin(q) + £ _1 cos(q) . 

Rearranging: 

tan(q) = 1/E. 

The above analysis assumes the sign of q is positive, which is intuitive for cruise flight. The 
singularity point where (sin (y) + £ _ cos(y)} = 0 defines the best glide angle (tan (y) = -1 /(L/D)) 
for T = 0, which in turn is derived from equation (A-5). 
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APPENDIX B 


OPTIMAL THRUST ANGLE FOR LANDING ROLLOUT 

The key parameter to be optimized in the landing rollout condition is the distance required to 
decelerate from touchdown to a stopped aircraft. This parameter can be expressed as: 

R = ^Vdt, (B-l) 

where the velocity, V, is expressed as 

V = V touchdown + j A x dt (B ' 2) 

and the deceleration, A x , developed in the “Optimal Thrust Angle and Benefits for Landing Rollout” 
section and expressed in equation (25), is: 

A x = g[rcos(ri)-D-^(W-rsin('n))]/W. (B-3) 

To find the maximum of A x with respect to r \ , the necessary condition is that the partial 

derivative dA /dr\ = 0, because maximizing the deceleration will minimize the landing rollout 
distance.^. Thus, dA x /dr\ = 0 = [g[(7’) - sin(r() - (0) - (0) + (p.r)(cos(r)))]]/W . Upon 
reducing, 0 = - sin(ri) + |X cos (r] ). Rearranging, tan(r|) = ji. 

The sign of T] is positive, which provides an additional down force on the wheels when thrust is 
reversed, thus increasing the rolling friction. The final result is the same for takeoff if the rotation 
velocity is assumed constant; however, as shown in this report, the vertical component of thrust reduces 
lift required for takeoff and thus the rotation velocity. 
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